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We consider the thermal Casimir effect in systems of parallel plates coupled to a mass-less free 
field theory via quadratic interaction terms which suppress (i) the field on the plates (ii) the gradient 
, of the field in the plane of the plates. These boundary interactions correspond to (i) the presence 

of an electrolyte in the plates and (ii) a uniform field of dipoles, in the plates, which are polarizable 
in the plane of the plates. These boundary interactions lead to Robin type boundary conditions 
£S) , in the case where there is no field outside the two plates. In the appropriate limit, in both cases 

Dirichlet boundary conditions are obtained but we show that in case (i) the Dirichlet limit breaks 
down at short inter-plate distances and in (ii) it breaks down at large distances. The behavior of 
the two plate system is also seen to be highly dependent on whether the system is open or closed. In 
addition we analyze the Casimir force on a third plate placed between two outer plates. The force 
acting on the central plate is shown to be highly sensitive to whether or not the fluctuating scalar 
field is present in the region exterior to the two confining plates. 
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I. INTRODUCTION 



The Casimir effect is often described in terms of how a boundary condition modifies the fluctuations of a field [l], [|J , 
the classic example being the case of the modification of the vacuum energy of the electromagnetic field between 
two conducting plates. However boundary conditions often arise from the consideration of ideal media such perfect 

S conductors. In reality the Casimir force is generated by interactions of the plates via the electromagnetic field, the 
material properties of the plates being coupled to the field. This point of view is embodied in the Lifshitz formulation 
of van der Waals interactions between macroscopic bodies @ . Also in the study of the critical Casimr force energetic 
boundary terms arise naturally in spin models due to surface interactions and fields 0, [E] . In this paper we analyze 
the fluctuation induced interactions due to a free mass-less field theory in the presence of planes with quadratic 
i interactions in the field variable. In an electrostatic analogy one case is equivalent to the field interacting with dipoles 

confined to the plane and the other case is equivalent to an electrolyte, in the Debye Huckel limit, confined to the 
<N| plates. In the limits where the dipole polarizability becomes infinite or the concentration of electrolyte becomes 
infinite, the limiting boundary conditions are Dirichlet. Clearly these two limiting cases mimic a conducting plate 
limit but via two distinct physical mechanisms. Here we show that the limit in which the Dirichlet limit is valid, 
for large but finite dipole or electrolyte strengths, depends on the model. In the electrolyte case deviations from the 
, Dirichlet limit become apparent at short inter plane separations but in the dipole case deviations appear for large inter 
plane separations. We also compare the results for two planes where the field exists in the region outside -the open 
system- with the case where the field does not exist outside - the closed system. This latter case corresponds to that 
arising in studies of the critical Casimir effect where the order parameter field exist within the critical fluid but not 
outside the boundaries of the physical system. The Casimir force in this case can be attractive or repulsive depending 
• • i on the boundary conditions at the two confining plates. However we show that when the fluctuating medium exists 
, ' outside the two plates then the interaction is always attractive. The case of a third plate confined between two other 
plates is also studied in both the open and closed systems. Here the force acting on the third plate can be evaluated 
' and we find a rich behavior and striking qualitative differences between the force on the central plate in the closed 
and open system. The method we use to carry out the computations is based on a path integral method adapted to 
planar geometries introduced in Q. The computations in this formalism are very short and straightforward and also 
have the advantage of highlighting immediately the differences between open and closed systems. 

II. THE MODEL AND THE TWO PLATE INTERACTION 

We consider a free scalar field theory analogous to that occurring for electrostatics with a free kinetic term every- 
where in space but with additional interaction terms with two surfaces at z — and z — I. 

H=\f dx[V0(x)] 2 + f dx6(z)M<l>(x)}+ I dx*(2-l)/ 3 [#x)]. (1) 
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The terms fi and f% are functionals of the field cf> on the two surfaces. Here we distinguish between the coordinate z 
perpendicular to the plates and the the coordinates perpendicular to the z direction denoted by x^. In this notation 
therefore any point is given by coordinate x = (z. xi). The field cj> is defined on a region of space with xj_ in a plane 
of area A and z in the region [— L, L] and we will be interested in the thermodynamic limits as A — > oo and L — > oo. 
Note that the open system corresponds to the so called defect plane case [4] , as opposed to the usual case considered 
in boundary critical phenomena where the field only exists in the region [0, 1] - physically the two cases are quite 
different as we shall see when comparing our results to some results in the literature 0, 0, HI, 0, [l0| . 

We will consider two types of interaction terms. First the case where the field 4> is acquires a mass in the plates 
(type I), i.e. it has a harmonic self-interaction 



/i[0(x)] = |^(x),. 



(2) 



This sort of interaction, for Cj positive, arises naturally in the Debye Hiickel theory of electrolytes and the coefficient 
Cj is proportional to the electrolyte concentration (for example see [6] ) When the Cj are positive this term will suppress 
the amplitude of the field <f> at the plates and we expect that in the limit Ci — > oo we will recover Dirichlet boundary 
conditions. We could also consider the case where the gradient of the field 4> i n the plane is energetically suppressed 
(type II) via 



/#(x)] = | [Vl^x)] 2 , 



(3) 



in this case Ei = — iVj_</> is suppressed and is set to zero in the limit Xi ~ > °°- This is the boundary condition for an 
electric field on a conductor. Clearly in both cases (up to an irrelevant zero mode) the boundary conditions for the 
two cases become equivalent in the limit Cj, Xi ~ > 00 ■ The purpose of this paper so to explore the modifications of the 
Casimir the effect when the coefficients Cj are finite. The boundary interaction term in Eq. ^ actually occurs quite 
naturally in a model of surfaces containing dipoles whose dipole moments are constrained to lie within the plane of 
the plates. The electrostatic Hamiltonian is now given by 

H = \ [ dx[V0(x)] 2 +i f dx6(z)V ±<l>(x) •Pi(x ± ) + i f dxS(z - i)Vj>(x) • P 2 (x_l) 
2 Jv Jv Jv 

+ -L / d x± -P 2 (x ± ) + -^ /dx ± -P 2 (x ± ) 
2xi J A 2x2 a J 



(4) 



where here —iV<j) is the electric field and Pi, 2 represent uniform dipole fields, constrained to lie within the plane of 
the plates, of polarizabilities Xi,2- Now integrating the corresponding partition function over the fields Pi and P2, 
yields an effective Hamiltonian for the field <f> with surface interaction terms of the form Eq. §3§ . 

We note that the classical equation of motion for the field in the electrolyte case induces Robin type boundary 
conditions at the surface relating the jump of the field derivative in the z direction to the value of the field on the 
surface. In the case where there is no space in the region external to the two plates, standard one sided Robin 
boundary conditions are obtained and this case has been extensively studied in the literature @, H, B E3] • 

In this paper we will use a calculational technique, based on the Feynman path integral method which has been 
adapted to study a variety of problems in electrostatics systems, the thermal Casimir effect and membrane fluctuations 
We proceed by decomposing the field (p into its Fourier components in the plane of x^, i.e. we write 



-7= X! 4>(K z) exp(ik • xi). 
V A , 



(5) 



The Hamiltonian is now given by 



dz 



d<j>Qs., z) 



k 2 0(k,z)0(-k,z) 



dz dz 

In the case of the scalar interaction (type I) term of Eq. ^ we have that 

5i(k) = Ci, 

and in the case of the transverse field (type II) interaction term of Eq. ([3]) we have 



i 5l (k)0(k,O)0(-k,O) + i 3 i(k)0(k,Z)0(-k,Z) 



(6) 

(7) 
(8) 
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where k — |k|. The resulting field theory is non interacting and the modes are all decoupled; we may thus write the 
partition function as a product over the partition function of the modes 



ln(Z) = ^ln(Z k ), 



(9) 



with 



Zk = J d[X k ]expf-^ J dz 



dXk' 
dz 



+ k 2 X 2 k 



9l (k)X k (0) 2 - L 2 (k)X k (lf 



(10) 



where we have decomposed the field <j> into its real and imaginary parts and as usual we only take half the sum over 
the modes k as the field <f> is real. Each partition function has the form of a simple harmonic oscillator path integral 
with interaction terms inserted at the times z = and z = I. The path integral kernel defined as 



pX(z 2 )=y 

K(x, y, zi, Z2, u, M) = I d[X]exp 



i)=x 



m r Z2 



dz 



dX' 

dz 



L0 2 X 2 



is given explicitly as 

K(x,y,z 1 ,z 2 ,u>,M) 



27rsinh(a>(z 2 — Z\)) 



exp ( — -(x 2 + y 2 )Mu> coth(uj(zi — z 2 )) + xyMLucosech(uj(zi — z 2 )) 



We now note that in the limit (z 2 — Z\) — > oo 



(11) 



(12) 



K(x,y,z 1 ,z 2 ,u,M) w I ——J 



1 



CX P(-77 W ( Z 2 - z l)) ex P 7r( X + V ) • 



Moj . 



(13) 



Thus the initial and final positions become decoupled. Therefore in the limit L — > oo, up to arbitrary terms depending 
on the values of the field x{—L) and x(L) we find 



Z k = exp(-kL+ -hi) J dxdy exp ^-^x 2 (k + gi{k))j K(x,y,0,l,k,(3)exp 



a 2 



y\k + g 2 {k)) 



(14) 



This is a trivial Gaussian integral to do and we find that, up to bulk terms denoted here by B k (independent of I and 
the gi), we have 



ln(Z fc ) =B k - X - 



ln(2fc + gi{k)) + ln(2fc + g 2 (k)) + In 1 



gi(fc)ff2(fc) exp(-2fc/) 
(2k + 9l (k))(2k + g 2 (k)) 



(15) 



The first term, as mentioned above, is a bulk term, the first term in the square bracket is a surface energy term 
for each surface and the final term is the I dependent term giving rise to the Casimir interaction. The I dependent 
Casimir free energy is thus given, in a space of total dimension d, by 



Fo(l) 



k B T 



(2tt)« 1 - 1 
k B T 



In 1 



gi (fc)g 2 (fc)exp(-2fcQ 
(2k + gi (k))(2k + g 2 {k)) 



(An)^T(^) 



k d - 2 dk\n{ 1- 



g\(k)g 2 (k) cxp(-2kl) 
(2k + 9l (k))(2k + g 2 (k)) J ' 



(16) 



where T is the Euler gamma function and the subscript o is to remind us that this is result for an open system. 

The first thing to notice is that in the strict limits g\ — > oo and g 2 — > oo we recover the classical result for Dirichlet 
boundary conditions 



F D (l) = k B TT(d-l)((d) 
A " (167r)^ 1 r(^i)/''- li 



(17) 
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where £ is the Riemann zeta function. In the limit where one of the gi is zero then the result is zero as it should be 
- this is a critical difference between the case of defect planes where the field exists outside the interior of the plates 
and the case where it does not exist outside the plates. We note that in a finite or closed system, where the field (f> 
does not exist outside the two plates as is the case in studies of the critical Casimir force @, H, H, > the result is 
somewhat different. Indeed in this case it is possible to have repulsive as well as attractive regimes, and moreover 
there is a residual interaction even in the case where one of the Ci is set to zero. In this open case the interaction is 
always attractive and it vanishes when either of the Ci, not just both, is set to zero. The results of @, H, @, El can 
easily be recovered in our formalism. When there is no region exterior to the slab the external propagators are absent 
and thus the ground state wave function at each interface is not there. This means that gi (which is added to k in our 
case) is simply replaced by gi — k. Upon subtraction of the bulk pressure one thus finds a Casimir force for a closed 
(hence a subscript c) system given by 

Fo(l) _ k B T f J d _ 2dk ^ / _ (, 9l (fc)-fc)(.g 2 (fc)-fc)exp(-2fc/) \ ^ 



A ( 47r )^ i r(^i) J V (k + 9l (k))(k + g 2 (k)) 

in agreement with the results of @, @, d, [l(| • Note that it is the appearance of the terms gi — k in the above expression 
that give the possibility of repulsive Casimir interactions 0, H, H U3| • The appearance of a repulsive interaction is 
most easily seen in the limit g\ — > oo and 172 — > 0. However we reemphasize that the presence of the field in the 
exterior region ensures that the interaction is always attractive (for gi positive). 

We now return to the case where the gi are finite, for the surface interaction term of Eq. ^ (type II) we find that 



Fq(1) _ ksT f d _ 2 ( _ X iX2fc 2 exp(-2fcQ 

A ( 4^^r(^)i V (2 + Xl k)(2 + X2 k) 



(19) 



Clearly in the large / limit the integral above is dominated by the small k behavior and thus for sufficiently large I 
the asymptotic behavior of the free energy is given by 

F Q)_ k B T f kd -, dkl J^X^eM-2kl)y (2Q) 



.4 



(4tt)— r(^) 



Thus at sufficiently large / the Dirichlet limit is no longer valid and the Casimir free energy will become dependent 
on the Xi ' I n this limit we thus find 

F il) _ k B TT(d +l) X iX2 (21) 



A i6(16 7 r) £ 2 i r(^)Z d + 1 ' 

and thus the strength of the interaction is considerably reduced. The large k part of the integral dominates the short 
distance behavior and the interaction thus remains of the Dirichlet form in this regime. The cross over length between 
Dirichlet and this modified long distance behavior is given by l c ~ X if the two X i are of the same order. If X i an d 
X2 are very different then there is an even richer behavior and it is possible to have an intermediate regime where 
F/A~-l/l d . 

Now we consider the case of the (type I) surface interaction Eq. ((2]), here we find that the small k limit agrees with 
the Dirichlet limit and thus the long distance behavior of the interaction in this case is of the Dirichlet form. The fact 
that the Dirichlet limit for type I interactions holds at large I is a consequence of the fact that Ci = 00 is a infrared 
stable fixed point (for both the free and interacting field theories) [H • However the deviations from the Dirichlet case 
are seen for large k and thus will show up in the short distance behavior of the interaction. In this case the Casimir 
pressure is given by given by 

p (l) = d F jl) = k B Tcic 2 r uddu exp(-u) 

° U dl A 2(1670^1X^)^-1 J iu + c 1 l){u + c 2 l)-c 1 C2l 2 exp(-u) [ ' 

In the limit of large I the Dirichlet limit is clearly always good, however it breaks down at small I when I <C 1 jc% . In 
this limit of small I we obtain (for d > 2) 

W 2(167r)^r(^±)^-i v ' 

It is easy to verify that this is a reduction of the Casimir pressure with respect to the ideal Dirichlet case in the 
limiting region where it is valid. 
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III. THE THREE PLATE INTERACTION 



In order to further demonstrate the power of the path integral method in the context of Casimir interaction we will 
consider the case of three plates. We keep two plates (plate (1) and (3)) at z = and z = I and we will place another 
plate between them at z = m. Again we denote the quadratic surface interaction coefficients by gi(k) where i is the 
plate number. The computation for this case within the path integral formalism is immediate (it encodes to a certain 
extent the transfer matrix formalism developed for van der Waals interactions in slab geometries developed in 13]). 
The partition function for the mode is given by 

Zk = (~J exp(-kL+ ^kl) J dxdydz exp ^-^x 2 (k + gi(k))J K(x, y, 0, to, k, (3) exp ^-^J/ 2 ff2(fc)) 

x K(y,z,m,l,k,l3)exp(-^z 2 (k + g 3 (k))) . (24) 



This yields 



\n(Z k ) = B k - X - [ln(2fc + gi (k)) + ln(2fc + g 2 {k)) + ln(2fc + g 3 (k))} 
1 ln[(l 9i{ k )92(k) exp(-2fcm) g 2 (k)g 3 (k) exp(-2fcm') 



2 LV (2k + 9l (k))(2k + g 2 {k)) (2k + g 2 (k))(2k + g 3 (k)) 
9i(k)gz{k)(g 2 {k) - 2k) exp(-2fc(TO + to')) 

(2k + gi(k))(2k + g 2 (k))(2k + g 3 (k)) U 1 ' 

where to' = I — m. The first term is a bulk energy independent of to and to', the second corresponds to three 
independent surface energies and the third contains the geometry dependent interaction. An important test of the 
above is that upon setting g 2 — we recover the two plate result of Eq. (fT5|) . The I (geometry) dependent part of 
the Casmir free energy is given by 

F (m,m') _ k B T f,d-2,,i m gi(k)g 2 (k) exp(-2km) g 2 (k)g 3 (k) e^q>(-2km') 



/k d ~ 2 dk lnlYl — 
i Ll (2k + gi (k))(2k + g 2 {k)) (2k + g 2 (k))(2k + g 3 (k)) 

9i{k)g 3 (k)(g 2 (k) - 2k) exp(-2fc(TO + to')) 

(26) 



(2k + gi(k))(2k + g 2 (k))(2k + g 3 (k)) 

If we take the Dirichlet limit gi — > for all i we obtain that the free energy is given by 

F D (m, to') = F D (m) + F D (m'), (27) 

i.e. the sum of the free energies of two independent systems with Dirichlet boundary conditions whose values are given 
by Eq. (|17p . This result is clearly expected on physical grounds as strict Dirichlet boundary conditions effectively 
decouple to two systems (plate 1 and 2 and plate 2 and 3). However in the general case we see that there is no 
decoupling and that n-body (plate) interactions are important. We also notice that the free energy also becomes 
equal to the sum of two independent terms (one dependent on to and the other on m') in the limit where g 2 — * oo. 

The case where the system is closed (no exterior field) can also be analyzed as before. Here we find (simply by 
replacing gi, 3 by gx >3 — k and leaving g 2 unchanged) 

F c (m,m') _ k B T f k d ~ 2 dkln[(l (di(k) - k)g 2 (k) esq>(-2km) 



(28) 



A (^) d -^T(^)J (k + 9l (k))(2k + g 2 (k)) 

g 2 (k)(g 3 (k) - fc)exp(-2fcTO') ( 9l (k) - k)(g 3 (k) - k)(g 2 (k) - 2k) exp(-2fc(TO + to')) 
(2k + g 2 (k))(k + g 3 (k)) (k + 9l (k))(2k + g 2 (k))(k + g 3 (k)) 1 



We see that as long as gi and #3 are finite then the results for the open and closed systems are quite different. 

Let us consider the case of type I boundary terms. For the case where the two outermost plates are fixed at a 
distance 1 let us define by 

V c , (x) = F Ct0 (x, l-x) = A—^——v c , (x), (29) 

(4tt) — T(^i) 




FIG. 1: Effective potential felt by a plane in the middle of two fixed planes all with type I boundary interactions. Solid lines 
for open systems and dashed lines for closed systems 



the effective potential felt by the central plate (plate 2). We restrict ourselves to the symmetric case c\ = C3 and 
which we will vary and we take c 2 = 1, also we shall consider the case d = 3. Shown in Figs. (la,b) are the scaled 
effective potentials v (evaluated by numerical integration) for the cases of open (solid line) and closed systems for 
Ci = 2 (a) ci = c 3 = 10 (close to the Dirichlet limit for the external plates). We see that for c\ = c 3 = 2 and 
c 2 = 10 that for an open system the middle plate is always attracted to the exterior plates. However for a closed 
system for c = 2 the middle plate is repelled from the two exterior plates and actually has an equilibrium position 
at the center of the two plates. For c — 10 the closed system has a potential which is close to that of the open 
system near the middle of the two plates and the midpoint is an unstable equilibrium point in both open and closes 
systems. However the closed system develops a repulsive potential close to the plates leading to a stable potential 
minima close to each plate. Notice that in the case of c\ — 10 that the deviations from Dirichlet behavior for the 
closed system are manifested when the distance between the central plate and the closest bounding plate is small, this 
should be expected from our discussion in section (II). If we consider the case of type II boundary terms we expect 
that deviations from the Dirichlet limit occur at large distances, We therefore consider a system of three plates again 
with a distance of 1 between the bounding plates. This distance should be large to see an effect and this is achieved 
by setting the polarizabilites x to be small. Shown in Figs. (2 a, b) are the scaled effective potentials v for the open 
(solid lines) and closed (dashed lines) systems. In Fig (2a) we have set xi — X3 — 0-1 and \2 = 0.1. In the case of 
an open system the central position is unstable and the middle plane is attracted towards the outer-plates. However 
for a closed system the central point is metastable and there is an energy barrier which must be crossed to reach 
the walls (which are ultimately attractive). If \i = (X3) is increased the local minima at the midpoint eventually 
disappears, as shown in Fig (2b) where we have taken xi = 0.3, and the two curves for the open and closed systems 
are qualitatively the same. 




IV. CONCLUSION 



In this paper we have studied a free field scalar theory in the presence of planes which interact quadratically with 
the field. In one case (type I) the field acquires a mass on the plane which suppresses its fluctuations. This would 
correspond to the way in which an electrolyte confined in the plane interacts with the thermal fluctuations of the 
electrostatic field. The second term (type II) is proportional to the square of the in plane gradient and arises due 
to dipole interactions with the electrostatic field. In both of these cases if the strength of the interaction is taken 
strictly to infinity we obtain Dirichlet boundary conditions. We have seen however that for finite interactions the 
interaction between the two plates deviates from the Dirichlet behavior, at small distances for type I and at large 
distances for type II and no longer has a universal form. We have also seen that for finite interaction terms there is 
a clear difference between open systems (where the fluctuating field exists outside the two plates) and closed systems 
(where there is no fluctuating field outside the plates). Notably for open systems the interactions between plates 
are always attractive, this is in contrast to the case of closed systems where it has been long established that both 
attractive and repulsive interactions are possible @, H, H, OH • We have also examined the behavior of a third plane 
sandwiched between two other planes, this demonstrates clearly the power of the path integral method used to to 
analyze Casimir-like interactions in planar systems. Again whether or not the system is open or closed can have a 
drastic influence on the force experienced by the third (central) plane. In closed systems the force felt by the central 
plane can be attractive or repulsive and even change sign in the same system, having positions of local equilibria away 
from the bounding walls (both stable and metastable). There are clearly many other configurations and set ups that 
one can study with the formalism developed here and it is possible that some of the basic mechanisms seen here can 
be exploited in the design of nano-devices [13] where Casimir type forces such as van der Waals interactions play an 
important role. 
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